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ABSTRACT 

We study the statistical properties of the combined emission of a population of discrete 
sources (e.g. X-ray emission of a galaxy due to its X-ray binaries population). Namely, 
we consider the dependence of their total luminosity Ltot = ^ Lk and of fractional 
rms tot of their variability on the number of sources n or, equivalently, on the nor- 
malization of the luminosity function. We show that due to small number statistics a 
regime exists, in which Ltot grows non-linearly with n, in an apparent contradiction 
with the seemingly obvious prediction (itot) — J L ^ dL cx n. In this non-linear 
regime, the rms tot decreases with n significantly more slowly than expected from the 
rms cx l/\/n averaging law. For example, for a power law luminosity function with a 
slope of Of = 3/2, in the non-linear regime, Ltot oc and the rms tot does not depend 
at all on the number of sources n. Only in the limit of n ^ oo do these quantities 
behave as intuitively expected, Ltot oc n and rms tot esc 1/ -^/n. We give exact solutions 
and derive convenient analytical approximations for Ltot and rms tot ■ 

Using the total X-ray luminosity of a galaxy due to its X-ray binary population 
as an example, we show that the Lx— SFR and Lx — relations predicted from the 
respective "universal" luminosity functions of high and low mass X-ray binaries are 
in a good agreement with observations. Although caused by small number statistics 
the non-linear regime in these examples extends as far as SFR < 4 — 5 Mq / yr and 
log(M*/Af0) < 10.0 - 10.5, respectively. 

Key words: methods: statistical - methods: data analysis - X-rays: galaxies - X-rays: 
binaries 



1 INTRODUCTION AND QUALITATIVE 
CONSIDERATION 

In many astrophysical situations a problem arises to pre- 
dict or interpret results of measurements of the total (com- 
bined) luminosity of a population of discrete sources. Among 
many examples are the total luminosity of X-ray binaries in 
a galaxy, or total flux of background sources detected above 
the sensitivity limit inside the fleld of view of a telescope. 

In the following discussion we will use high mass X-ray 
binaries (HM XB) in star forming galaxies as an ex ample. As 
was shown bv lGrimm. Gilfanov fc Sunvaev^^l2003^ . the lumi- 
nosity distribution of HMXB sources in a galaxy is described 
to first approximation by a "universal" luminosity function 
whose shape is the same in all galaxies and whose normal- 
ization is proportional to the star formation rate (SFR) of 
the parent galaxy: 

^ cx SFR X /(L) (1) 
In a broad luminosity range, log(Lx) ~ 35.5 — 40.5, the 



shape of the HMXB "universal" luminosity function is close 
to a power law, f{L) = L~°', with the slope of a « 1.6. 
Importantly, the luminosity of the compact sources in star 
forming galaxies appears to be restricted by a maximum 
value of Lcut ~ few x 10*° erg/s. This cut-off luminosity 
can be defined, for example, by the Eddington luminosity 
limit for the most massive objects associated with the star 
forming regions. Obviously, on the faint side, the luminosity 
distribution eq.Q must become fiatter or have a cut-off as 
well, in order to keep the total number of sources finite. 
This low luminosit y cut-off may be caused for example by 
thepropeller effect ('lUarionov fc SunvaevllOTsI) as discussed 
byihtvkovskii fc Gilfanov ( 2004) . 

The expectation value for the total number of sources 
in a galaxy equals 

(iVtot) = ^ — dL cx SFR (2) 

and, naturally, is directly proportional to its star formation 
rate. The number of sources actually observed in a given 
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galaxy obeys Poisson distribution p^{N) with n defined by 
eq.l|5J. Apart from effects of counting statistics, tfie number 
of HMXB sources found in an arbitrarily chosen galaxy will 
be close to the above expectation value. 

The problem, considered in this paper, is the behavior of 
the total luminosity of high mass X-ray binaries in a galaxy 



k = N 

Ltot = ^ ] Lk 
fc=i 



(3) 



as a function of its star formation rate. 

An apparently obvious expression for the Ltot can be 
obtained integrating the luminosity distribution Q: 



(Ltot) 



dN 

L — dL oc SFR 
dL 



(4) 



Hence, one might expect that the total X-ray luminosity 
of HMXB sources is also directly proportional to the star 
formation rate of the host galaxy, as is the total number 
of sources. This problem, however, involves some subtleties 
related to the statistical properties of the power law dis- 
tribution of the sources over luminosity, which appear not 
to have been recognized previously, at least in astrophysi- 
cal cont ext (a somewhat rela ted problem has been consid- 
ered by iKaloeera et"al]l200ll in connection with estimating 
the coalescence rate for the NS-NS binaries in the Galaxy). 
Although eq.Q correctly predicts the average X-ray lumt- 
nosity computed for a large number of galaxies with similar 
values of star formation rate, it fails to describe the relation 
between the most probable value of X-ray luminosity of an 
arbitrarily chosen galaxy and its SFR. The main surprise of 
the study presented here is that, in the low SFR regime, the 
relation between SFR of the host galaxy and the total lumi- 
nosity of its HMXBs is non-linear - with increase of the star 
formation rate the luminosity appears to grow faster than 
linear. The relation becomes linear only for sufficiently high 
star formation rates, when the total number of objects with 
a luminosity close to the maximum possible value, defined 
by I/cut, becomes sufficiently large. 

This can be illustrated by the following simple consider- 
ation. For an arbitrarily chosen galaxy, the brightest source 
is most likely going to have a luminosity Lmax defined by 
the condition^ 



dN 



(5) 



For a power law luminosity distribution with slope a and 
with a cut-off at Lcut, eq.Q, the above expression yields: 



I/max oc SFR — low SFR 
Lmax = Lcut high SFR 



(6) 



As might be intuitively expected, at low SFR, the 
most probable luminosity of the brightest source increases 
with SFR, until it reaches the maximum value of Lcut. 
The threshold value of the star formation rate, separat- 
ing low and high SFR regimes, is defined by the condition 
N{L ~ Lcut) ~ 1, i.e. that there are ~few sources expected, 
with luminosity close to the cut-off value Lcut. 

^ Indeed, for example sources 10 times more luminous will appear 
on average in one out of ~ 10'"~^' galaxies. 



The most probable value of the total luminosity, Ltot 
can be then computed integrating the luminosity function 

from Lmin to Lniax^ 



dN 

11' 



L dL 



which, for 1 < a < 2 and Lmin << Lmax, leads to 
SFR^ low SFR 



Ltot oc 



SFR high SFR 



(7) 



(8) 



i.e. is non-linear in the low SFR regime and becomes linear 
only at high star formation rates. 

This can be qualitatively understood as follows. For a 
slope of the luminosity distribution 1 < a < 2, the total 
luminosity of a galaxy is defined by the brightest sources. 
The non-linear behavior in the low-SFR limit is caused by 
the fact that an increase of the SFR leads to non-linear in- 
crease of the luminosity of the brightest sources. Therefore 
their total luminosity grows faster than the star formation 
rate. This non-linear growth continues until the maximum 
possible value of the luminosity of the compact sources is 
achieved. Further increase of the star formation rate leads 
to a linear increase of the number of the brightest sources 
in the galaxy, but not of their individual luminosities. Con- 
sequently, the Ljf — SFR relation becomes linear. 

In the more formal language of statistics such a behav- 
ior is related to the properties of the probability distribution 
of the collective luminosity p(Ltot). In particular, it can be 
understood in terms of the difference between the expecta- 
tion mean and the mode of the probability distribution. The 
expectation mean is defined as 



(Ltot) 



p (Ltot) Ltot dLt, 



(9) 



and is given by eq.Q. The mode of the statistical distribu- 
tion, Ltot , is defined as the value of the random variable (Ltot 
in our case) at which the probability distribution p(Ltot) 
reaches its maximum value. Whereas the expectation mean 
(Ltot) describes the result of averaging of the X-ray lumi- 
nosities of many galaxies having similar values of SFR, it is 
the mode of the p(Ltot) distribution that predicts the most 
probable value of the total luminosity of a randomly cho- 
sen galaxy. The non-linear behavior in the low SFR regime 
is caused by the skewness of the probability distribution 
p(Ltot) resulting in a difference between expectation mean 
and mode. In the high SFR limit, the p(Ltot) distribution 
asymptotically approaches the Gaussian distribution, in ac- 
cord with the Central Limit Theorem. The boundary value 
of SFR, separating the non-linear and linear regimes of the 
Lx— SFR relation is defined by the parameters of the lumi- 
nosity function. 

Interestingly, the fact of the existence of a linear regime 
in the Lx-SFR relation is a direct consequence of the cut-off 
in the luminosity function. Only in the presence of a maxi- 
mum possible luminosity for the sources, Lcut, (for instance 
Eddington limit for a neutron star) a linear regime can ex- 
ist, when the total luminosity of a galaxy is defined by a 
sufficiently large number of bright sources near Lcut, and 
any subsequent increase of the star formation rate results in 
linear growth of the total luminosity. 
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In the above discussion, we used high mass X-ray bi- 
nary populations in star forming galaxies as an exam- 
ple. Obviously, the effect considered in this paper is of a 
broader general interest and is at work in many situations 
related to computing/measuring integrated luminosity of 
a finite number of discrete sources. One example related 
to the application of stellar synthesis models to observa- 
tio ns of stellar clus ters has been independently discovered 
bv lCervino fc Luri dana ( 200J). 

The paper is structured as follows. In section|5|we con- 
sider the statistical properties of the total luminosity, in par- 
ticular in section 12.11 derive the formulae for the probabil- 
ity distribution p (Ltot), using two different approaches, and 
present results of numerical calculations in 12.21 The vari- 
ability of the total emission is studied in section |3] In sec- 
tion |1| we discuss astrophysical applications, including the 
properties of the total X-ray emission of a galaxy due to its 
population of high- (section 14.21 14.311 and low- (section 14.411 
mass X-ray binaries. Our results are summarized in section 
13 In the Appendixes we derive convenient approximations 
for Ltot and fractional rms of the total emission and con- 
sider their asymptotical behaviour. A casual reader, not in- 
terested in the mathematical aspects of the problem, can 
skip section inl and proceed with section 



2 TOTAL LUMINOSITY 

We consider a population of sources with a luminosity func- 
tion (LF): 



(10) 



The expectation values for the total number of sources and 
the total luminosity are: 



(n) 



dN 

H 



dL 



(11) 



of sources in the interval {Lk, Lu + 5Lk) obeys a Poisson 
distribution with mean 

tik = ^{L = Lk)5Lk (14) 

For 5Lk 0, it is sufficient to consider the occurrence of 
either zero or one source per interval whose probabilities are, 
respectively: 

dN 

11 

dN 
'dL 

The probability distribution for the combined luminosity of 
the sources in the fc-th interval is 

dN 



Pk{0) = l-^SLk+O{SLl] 



Pfe(l) = ^SLk + 0{SLl) 



(15) 



P (Ltot.k) = (5(Ltot,k) (^1 - -jj- SLk^ + 



(16) 



where 5{x — xq) is the delta function. The characteristic 
function of this distribution is 



(17) 



p(itot,k)= J p(I/tot,k) e*"-^ dL = 

dN ^dN 
1 - -77- 5Lk + -rr SLk e 
dL dL 

Using the convolution theorem, the characteristic function 
of the probability distribution of the total luminosity can 
be computed as a product of characteristic functions of 
p (Ltot.k) 



K-^-tot) = JJp(Ltot,k) 



lnp(Ltot)=^ln(l-^5L 



(18) 



dN ,r 

I- H oLh e ' 

dL dL 



poo , J. J 

{Ltot) = j L — dL (12) 
It is assumed that both quantities are well defined and finite. Finally 



2.1 The probability distribution of the total 
luminosity 

Below we present two methods of computing the probability 
distribution of the total luminosity. The second of these two 
methods is somewhat more convenient computationally. 



(itot) = / p(Ltot)e-'"^*°'d^: 



exp ( I e^^ dL — (n) — icuLtot ] du! 



where (n) is given by eg. lllll . 



(20) 



2. LI Method I 

To compute the probability distribution for the total lumi- 
nosity, p (Ltot), we divide the Li — L2 luminosity range into 
intervals of infinitesimal width SLk and express the com- 
bined luminosity of all sources as a sum 

Ltot = ^ ^ Ltot.k (13) 

where Ltot.k is the combined luminosity of the sources in 
the fc-th interval, running from Lfc to Lk + SLk ■ The number 



2 A. 2 Method II 

The probability distribution of the total number of sources 
follows a Poisson distribution P^{n) = e^''/n! with mean 
fi given by ea. lllll . The probability distribution of the total 
luminosity is: 

n— 00 

p(Ltot) = PAn)p4Ltot) (21) 

where Pn(Ltot) is the probability distribution of the total 
luminosity of exactly n sources. In the majority of practically 
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Figure 1. Probability distributions of the average luminosity, Lj^^/n, of n discrete sources with a power law LF, eq.J^. The lower 

and upper luminosity cut-offs were fixed at Li = 1 and L2 = 10^ for all plots. The value of the slope a is indicated in each panel. Each 
curve is marked according to the number of sources n. The vertical dashed lines show the expectation mean (Ltot) /n. 



interesting cases, the total number of sources is sufficiently 
large, (n) >> 1, and the Poisson distribution in ea. l|21|l can 
be replaced by the delta- function 5{n — (n)): 

p{Ltot) ~ Pn{Ltot), n^{n)»l (22) 

where n = (n) is given by eg. lllH 

The probability distribution of the sum of n random 
variables, Ltot = Sfe-i -^fe' '^^^ ^® calculated using the con- 
volution theorem: 



Pn = Pi 

poo 

pi{uj) = / pi{L)e^'^dL 
Jo 



(23) 



where pi (L) is the probability distribution for the luminosity 
of one source, which equals the luminosity function with 
appropriate normalization: 

1 diV 



The probability distribution Pn(-Ltot), defined by 
ea. l|28^ . describes the case when n sources are observed. On 
the contrary, the distribution p(Ltot), defined by eg. 12011 . is 
parametrized via the normalization of the luminosity func- 
tion or, equivalently, via the expectation value (n), and de- 
scribes the case when (n) sources are expected. These two 
distributions are related via ea. l2H . For n >> 1, which is 
often the case, they are nearly identical, and it can be as- 
sumed that n = (n) and both quantities are related to the 
normalization of the luminosity function via ea. dllll . 



2.2 Illustrative examples 

To illustrate the behaviour of the total luminosity we assume 
a power law luminosity function: 



dN 



AL-°' if Li < L < La 
otherwise 



(25) 
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Figure 2. The ratio of the most probable value of the total luminosity Ltot to the expectation mean (Ltot) versus number of sources 
observed in Li — L2 luminosity range for different values of the LF slope a and the ratio L2/L1. The results of exact calculation using 
ea. l2iji are shown by the solid symbols connected with thick line. The thin solid line shows the approximate relation calculated from 
eas. lA7i and iA3l . as described in the Appendix lAll 
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Figure 3. The most probable value of the total luminosity Ltot and its intrinsic dispersion versus number of observed sources for different 
values of the LF slope a. The lower and upper cut-off luminosities were fixed at Li = 1, L2 = 10^. The results of exact calculation 
using eg. 121^1 . Behavior of Ltot is shown by the thick solid line. The narrower and broader dashed areas correspond to its 67% and 90% 
intrinsic dispersion. The thin solid line shows linear relation for the expectation mean (Ltot) oc n. 



2.2.1 Probability distribution p„ (Ltot) 

The probability distributions p„{Ltot) for various values of 
the slope a and the number of sources n are shown in Fig0 
To facilitate the comparison of distributions for different val- 
ues of n, the abscissa in these plots is the average luminosity 
Ltot/n = X/fc=i Lk/n. The values of the luminosity cut-offs 
were fixed at Li = 1 and L2 = 10'^ 

The skewness of the probability distribution Pn{Ltot) 
leads to a deviation of its mode Ltot (the most probable 
value of the total luminosity) from the expectation mean 



{Ltot} indicated in each panel by the vertical dashed line. 
The effect is most pronounced for a ~ 3/2, is unimportant 
for for shallow luminosity functions with a < 1 and gradu- 
ally diminishes with increasing a at a > 2. 

For illustration, we also show in FigQthe probability 
distributions for a fiat LF, a = 0, in which case Ltot exactly 
equals (Ltot). Naturally, for any value of a the p„(Ltot) — > 
Gaussian distribution in the limit of n — > 00, in accord with 
the Central Limit Theorem. Correspondingly Ltot ^ (^tot) 
in this limit. 
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2.2.2 The most probable value of total luminosity 

The dependence of the most probable value of the total lu- 
minosity I/tot on the total number of sources n for different 
values of a and the ratio L2 / Li is shown in Fig|5|and|3 Con- 
venient analytical approximations for the Ltot ~ n relation 
are derived in Appendix^and its asymptotical behavior is 
considered in Appendix I A2I 

The Z/tot — n relation is significantly non-linear for 
"small" number of sources or, equivalently, for small values 
of the LF normalization A. For Li << L2, a > 1, the bound- 
ary between the non-linear and linear regimes, expressed in 
terms of the normalization A, depends only on the LF slope 
a and its high luminosity cut-off L2 (see Appendix IA2II . 
This is due to the fact that the behavior of Ltot is defined 
by the number of sources near the high luminosity cut-off 
of the luminosity function, rather than by the total number 
of sources in the entire Li — L2 luminosity range, i.e. the 
condition for the linear regime is: 



N{L ~ La) 



dN 



L2) >1 



(26) 



The total number of sources in the L1 — L2 luminosity range, 
on the contrary, is defined by the low luminosity cut-off Li 
(for a > 1). For sufficiently large L2/L1, the non-linear 
regime can occur for an arbitrarily large total number of 
sources (cf. the curves corresponding to different values of 
L2/L1 in Fig|2J. Interestingly, for the LF slope in the range 
of 1 < Q < 2, where the effect is strongest, there is a rel- 
atively sharp break separating the non-linear part of the 
dependence from the linear part. 

Although for small n the most probable value of lumi- 
nosity Ltot can deviate significantly from the expectation 
mean (Ltot) ('Fia lll2L the condition J p (Ltot) Ltot dLtot = 
(Ltot) is satisfied for any n. Consequently, the average of Ltot 
over a large number of realizations with the same n always 
equals (Ltot). This equality is achieved due to the existence 
of outliers, having values of Ltot significantly exceeding both 
Ltot and (Ltot), in accordance with the skewness of the prob- 
ability distribution pn(Ltot) for small n. This naturally leads 
to enhanced and asymmetric dispersion of the observed val- 
ues of Ltot in the non-linear regime as illustrated by the 
shaded areas in Fig|3| 



2.3 The luminosity of the brightest source in a 
sample 

The probability distribution for the luminosity of the bright- 
est source observed in a sample of n sources equals: 



p(Lmax) = [pi(L < Lmax)]" ^Pl(Lmax) n 



(27) 



where pi (L) is the probability distribution for the luminosity 
of one source (e.g. eg. 1241 ') and pi(L < Lmax) denotes the 
cumulative probability pi(L < Lmax) = /p'"'"'"' Pi(i) dL. 

The p(Lniax) distribution for a = 1.5 is shown in FigU] 
and illustrates the intuitively obvious fact that if the num- 
ber of sources is sufficiently small the brightest sources most 
likely will not reach the highest possible value of L2 . An ana- 
lytical formula for the most probable value of the luminosity 
of the brightest source is given by eq. lAH . 
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Figure 4. Probability distributions p(Lmax) of the luminosity of 
the brightest source in a sample of n sources with a power law 
LF with the slope a = 1.5 and Li = 1, L2 = 10^. Each curve is 
marked according to the number of sources n. 



3 VARIABILITY OF THE TOTAL EMISSION 

For a population of n sources with luminosities Lk and frac- 
tional rms of aperiodic variability rmsk, the fractional rms 
of the total emission equals: 



En J- 2 2 



(28) 



assuming that variations of the source fiuxes are uncorre- 
lated with each other. In the following we also assume for 
simplicity that all sources have the same value of fractional 
rms, i.e. rmsfc = rmso. 

In the limit of n 00, corresponding to the linear 
regime in the Ltot ~ n relation, the sums in the eg. (1281 can 
be replaced by the respective integrals of the LF: 



it _ jL^dN/dLdL 
4 (/ LdN/dLdL) 



1 



oc — oc — oc 
2 n A 



1 

Ltot 



(29) 



In the linear regime the fractional rms of the collective emis- 
sion obeys oc averaging law, as expected for uncorre- 
lated variations of individual sources. 

In the non-linear regime, however, for a sufficiently flat 
luminosity function, the total luminosity is defined by a few 
brightest sources. To first approximation the number of such 
sources effectively contributing to the sums in eg. 1281 does 
not depend on the LF normalization. Consequently, the frac- 
tional rms of the total emission is constant, independently 
on the total number of sources or of their total luminosity. 
This contradicts the intuitive expectation that the fractional 
rms decreases with the number of sources as rms oc 

To illustrate this behavior we performed a series of the 
Monte-Carlo simulations for a power law LF. For each set 
of parameters and a given value of the number of sources, 
in each run n sources were placed between Li and L2 
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Figure 5. Variability of the total emission. The ratio of the fractional rms of the combined emission to that of one source, rmso, as a 
function of the number of sources for different values of the LF slope a (indicated in the upper-right corner of each plot) and different 
values of the ratio of L2/L1. The thick solid lines show the square root of the most probable value of the rms'^^^^/rmsQ, the shaded 
areas show the 67% intrinsic dispersion, both obtained from the Monte-Carlo simulations. The thin solid lines were computed using the 
approximation given by eg. iBSi . The dashed line shows rms tot / T'msQ oc l/y/n dependence. 



with a power law luminosity distribution, eg. (1251 . For each 
run, the ratio rms\ot/rmsQ was computed following eg. 12811 . 
From results of many runs, the probability distribution for 
rms tot /'''tisq was obtained. The maximum of this distri- 
bution defines the most probable value of rmstot and its 
width characterizes the intrinsic dispersion of this quantity. 
The examples are shown in Fig|5]for two values of the LF 
slope, a = 1.5, 2.5. For q < 1 and a > 3 there is no non- 
linear regime and the fractional rms of the total emission 
obeys oc l/\/n law for any n. In plotting Fig^l we converted 
rms tot to rms tot, so that the thick solid curves shows the 
behavior of the square root of the most probable value of 
the rms tot /rmso ratio. 

Approximate formulae for rms tot are derived in Ap- 
pendix^ 



4 ASTROPHYSICAL APPLICATIONS 

4.1 Determining LF parameters from total 
emission of unresolved sources 

The shape of the Ltot — A relation {A is the normalization of 
the luminosity function) is defined by the parameters of the 
luminosity function. For 1 < a < 2, it has two distinct power 
law regimes, separated by a break (Figl^andl^: 



Ltot oc 



A 



A> Ay, 



(30) 



The position of the break between the non-linear and linear 
regime is defined by the high luminosity cut-off of the LF 
and its slope as described by eg. lA13j . This opens a pos- 
sibility to determine the LF parameters without actually 
resolving the individual sources, but studying large samples 
of objects (e.g. galaxies) and the relation between their total 



emission and the normalization A. The value of A in many 
cases can be determined independently from observations at 
other wavelengths. For example, the normalizations of the 
luminosity functions of high and low mass X-ray binaries 
are p roportional to the star formation rate iGrimm et alJ 
12003) and stellar mass (lGilfanovll200l of the host galaxy 
respectively. Both quantities can be determined from the 
conventional indicators, such as radio or near-infrared lumi- 
nosities. 

Of course with the sub-arcsec angular resolution of 
Chandra the luminosity distribution of point sources in 
nearby galaxies can be studied directly. However, for more 
distant galaxies, D > 30 — 50 Mpc, even Chandra resolution 
becomes insufficient and only the total luminosity of the 
galaxy can be measured. Provided that the contaminating 
contribution of the emission of a central AGN and/or of hot 
X-ray emitting gas can be constrained or separated, one can 
study the relation between the total luminosity Lx of a (un- 
resolved) galaxy and its star formation rate or stellar mass. 
The Lx— SFR or Lx — M« "growth curves" constructed for 
large samples of galaxies and spanning a broad range of SFR 
and M* can be used to constrain the XLF parameters of X- 
ray binaries in distant galaxies including those located at 
intermediate and high redshifts. With this, one can study 
the influence of a number of factors such as effects of binary 
evolution, metallicity, regimes of star formation, etc. on the 
luminosity distribution of X-ray binaries. 



4.2 High mass X-ray binaries in star forming 
galaxies 

4-2.1 Lx^SFR relation for star forming galaxies 

The slope of the "universal" luminosity function of IfMXBs, 
a = 1.6, is in the range where the non-linear behavior of the 
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Figure 6. Left: The probability distributions p(Ltot/SFR) for different values of SFR. The vertical dashed line shows the expectation 
mean, defined by eq.gj. Right: The Lx-SFR relation. The open circles are nearby galaxies observed by Chandra, the filled triangles 
are spatially unresolved nearby galaxies observed by ASCA and BeppoSAX, for which only total luminosity is available, the filled circles 
are distant star forming galaxies from the HDF-N. The thick grey line is relation between the SFR and the most probable value of the 
total luminosity, predicted from the "universal" HMXB XLF, the shaded area shows its 67% intrinsic spread, the dashed line is the 
expectation mean, defined by eq.j^. The five nearby galaxies, used bv iGrimm ct al. ( 2n03.'l to derive the "universal" HMXB XLF, are 
marked as crossed boxes. 



itot — A relation is most pronounced. In the left panel of 
FiglHl we plot the probability distribution of the total lu- 
minosity p(Ltot), computed for different values of the star 
formation rate. The distribution is strongly asymmetric in 
the non-linear low-SFR regime which, for the parameters of 
the "universal" HMXB XLF, corresponds to the formation 
rate of massive stars below SFR ^4 — 5 /yr. Note that 
a small value of SFR does not necessarily imply a small to- 
tal number of sources, which is defined by the (unknown) 
low luminosity cut-off of the HMXB XLF. For example, for 
SFR=0.2 M0/yr, when the non-linear effect is very pro- 
nounced, the total number of sources might be as large as 
~ 300 (~ 1200) for the low luminosity cut-off of lO^"* (10^^) 
erg/s. These low values of the star formation rate corre- 
spond to the familiar examples of the Milky Way galaxy 
and the Magellanic Clouds. On the opposite end among the 
relatively nearby and well known galaxies are the Antennae 
interacting galaxies which, with the star formation rate of 
SFR~ 7 M0/yr, have a nearly symmetric p(I/tot) distribu- 
tion, sufficiently close to the normal distribution. 

The predicted Lx— SFR relation is shown in the right 
panel in Fig|S| along with the measured values of X-ray 
luminosities and the star formation rates for a number of 
nearby galaxies and galaxies observed with Chandra at in- 
termediate redshifts, z ~ 0.2 — L3, in the Hubble Deep 
Field North. The data shown in Fig|S|are from lCrimm et all 
(I2OO3I), complemented with the local galaxies data from 
liianalh, Comastri & Seti (2003). In plotting the latter we 
removed the duplications and the galaxies likely to be con- 
taminated by the contribution of low mass X-ray binaries, 
unrelated to the current star formation activity, as dis- 
cussed bv lCilfanov et all ll2004f) . The luminosities and star 



formation rates for the Hubble Deep Field North galaxies 
jBrandt et all200ir) were computed for the following cosmo- 
logical parameters: Ho = 70 km/s/Mpc, Q.m = 0.3, A = 0.7, 
as described in Gilfanov et al. (2004). 

Fig|S|further illustrates the difference between the mode 
and the expectation mean of the p(Ltot) probability distri- 
bution. The thick solid line in the right panel shows the 
SFR-dependence of the mode of p(Ltot) and predicts the 
most probable value of the X-ray luminosity of a randomly 
chosen galaxy. If observations of many (different) galaxies 
with similar values of SFR are performed, the obtained val- 
ues of Ltot will obey the probability distribution depicted in 
the left panel. The average of the measured values of Ltot 
will be equal to the expectation mean given by eq.Q and 
shown by the dashed straight lines in the left and right pan- 
els. Due to the properties of p (Ltot) these two quantities are 
not identical in the low-SFR limit when the total luminosity 
is defined by the small number of the most luminous sources. 
Only in the large SFR limit when there are sufficiently many 
sources with luminosities near the cut-off of the luminosity 
function, log(Lx) ~ 40, does the Ljf — SFR relation behave 
in the intuitively expected way. 

Owing to the skewness of the probability distribution 
p(Ltot) large and asymmetric dispersion among the mea- 
sured values of Ltot is expected in the non-linear low-SFR 
regime. This asymmetry is already seen in Fig|S|- at low 
SFR values there are more points above the thick solid curve 
than below. Moreover, the galaxies lying significantly above 
the solid and dashed curves in Fig|S| should be expected at 
low SFR and will inevitably appear as the plot is populated 
with more objects. Such behavior differs from a typical astro- 
physical situation and should not be ignored when analyzing 
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Figure 7. Dependence of the Lx — SFR relation on the XLF cut- 
off luminosity Lcut- The three curves, corresponding to different 
values of Lcut coincide in the non-linear low SFR regime but differ 
in the position of the break between linear and non-linear regimes. 
The data are the same as in Fig|^ 

and fitting the Lx-SFR relation in the low SFR regime. In 
particular, due to non-Gaussianity of the p(Z/tot) distribu- 
tion the standard data analysis techniques - least square 
and fitting - become inadequate. 



Figure 8. Variability of the total emission of HMXBs in star 
forming galaxies. Dependence of the ratio rms tot /rmso on the 
star formation rate. The thick solid line shows the most probable 
value of the rms tot /rmsQ , the shaded area shows its 67% intrinsic 
dispersion, both obtained from the Monte-Carlo simulations for a 
power law LF with the slope of o = 1.6 and a cut-off at L2 = 2 ■ 
10*" erg/s. The dashed line shows asymptotical behavior at large 
SFR. The filled circles correspond to HMXB sources in the Milky 
Way, SMC and the Antennae galaxies, computed from eg. 1281 . 
using the observed luminosities of X-ray sources in these galaxies. 



4.2.2 High luminosity cut-off in the HMXB XLF 

The position of the break between the non-linear and lin- 
ear parts of the Lx-SFR relation depends on the LF slope 
and its cut-ofTIuminosity (Fig|7| eg IA13t : SFRbroak oc L"~^^. 
This allows to constrain the parameters of the luminosity 
distribution of compact sources using the data of spatially 
unresolved galaxies as discussed in section im 

Agreement of the predicted Lx— SFR relation with 
the data both in high- and low-SFR regimes confirms the 
universality of the HMXB luminosity function, derived by 
iGrimm et alJ i2003h from significantly fewer galaxies (shown 
as crossed boxes) than plotted in Figs [HI Q It provides an 
independent confirmation of the existence of a cut-off in 
the HMXB XLF at log(Lcut) ~ 40.5, including HMXBs in 
spatially unresolved high redshift galaxies from the Hubble 
Deep Field North. This implies, in particular, that ultra- 
luminous X-ray sources (ULX) at redshifts of 2 ~ 0.2 — L3 
were not significantly more luminous than those observed in 
the nearby galaxies. 

4.2.3 Aperiodic variability 

X-ray flux from X-ray binaries is known to be variable in a 
broad range of time scales, from ~msec to ~ years. In addi- 
tion to a number of coherent phenomena and quasi-periodic 
oscillations, significant continuum aperiodic variability is of- 
ten observed. The fractional rms of aperiodic variations de- 



pends on the nature of the binary system and the spectral 
state of the X-ray source and is usually in the range from a 
fraction of a per cent to ~ 20 — 30 per cent. Correspondingly, 
the combined emission of X-ray binaries in a galaxy should 
be also variable in a si milarly broad range of time scales. It 
has been suggested by iGrimm. Gilfanov fc SunvaevI (|2004) 
that due to a large difference in the characteristic time scales 
of the accretion flow onto a stellar mass object and onto a su- 
permassive black hole variability of the X-ray emission from 
a galaxy can be used to distinguish between the combined 
emission of a population of X-ray binaries and that of an 
accreting supermassive black hole in the centre of a galaxy 
(AGN). 

From results of section |21 and Appendix ^ one should 
expect that in the non-linear low-SFR regime the fractional 
rms of the total X-ray emission of a star forming galaxy is 
independent of SFR. This prediction is confirmed by the 
results of the Monte-Carlo simulations performed as de- 
scribed in section 13 and shown in Fig|H| For moderate star 
formation rates, SFR < 5 — 10 Mq /yr, we predict a rather 
large aperiodic variability of the total emission of HMXBs 
at the level of ~ 1/3 — 1/2 of the fractional rms of individ- 
ual X-ray binaries. At larger values of SFR, corresponding 
to the linear regime in the Lx— SFR relation, it decreases 
as rms oc l/vSFR, in accord with eg. lBSII . Also shown 
in FiglHl are values of the fractional rms reduction, com- 
puted directly from eg. 12811 using th e luminosities of th e 
observed HMXBs in the Milky Way JCrimm et alJl2002f) . 
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Figure 9. Left: The probability distribution of the luminosity of the brightest HMXB source for different values of the star formation 
rate, computed from eg. 1271 using the parameters of the "universal" luminosity function of HMXBs. Right: Expected luminosity of the 
brightest HMXB vs. star formation rate. The thick solid line shows the most probable value of the I/max feg. lAli '). the shaded area 
shows its 67% intrinsic dispersion, obtained from the probability distribution given by eg. 1271 . The filled circles show maximum observed 
luminosity of HMXBs in the Milky Way and several nearby star forming galaxies. 



SMC J_Yokogawa ct al."2000'), and in the Antennae galax- 
ies (iZezas et al 2002*1 . Note, that the predicted rms— SFR 

relation can be modified by the luminosity dependence of 
the rms of individual sources. This factor might become es- 
pecially important at large values of SFR when the total 
luminosity of a star forming galaxy is dominated by ULXs 
whose variability properties we know little about. 

4-2.4 Luminosity of the brightest source 

As the first Chandra observations of compact sources in 
nearby galaxies bec a me available, it has been noted (e.g. 
I Sarazin et all I 20011: llrwin et al.l 120021: iFabbiano fc White! 
l2003 h that the luminosity of the brightest X-ray binary in 
a galaxy might depend on its properties. In particular in 
the case of high mass X-ray binaries it appeared to corre- 
late with the star formation rate of the host galaxy. For 
example, in the Antennae galaxies, a number of compact 
source s have been disco vered with luminosities of ^ lO*"' 
erg/s (IZezaseTaOEool). On the other hand, the luminosi- 
ties of the brightest HM XB sources in the M ilky Way do 
not exceed < 10^* erg/s llCrimm et al.ll2002^ . It has been 
argued that this might reflect the difference in the intrinsic 
source properties, related to the difference in the galactic 
environment and in initial conditions for X-ray binary for- 
mation in starburst galaxies and in galaxies with weak and 
steady star formation. 

However, as discussed in section 12.31 the probability 
distribution for the luminosity of the brightest source in a 
galaxy, p(Linax), depends on the LF normalization, i.e. on 
the SFR of the host galaxy in the case of HMXBs. The lu- 
minosity of the brightest source increases with SFR, until it 



reaches the maximum possible value, deflned by the high lu- 
minosity cut-off of the LF. The p (I/max) distributions, com- 
puted from eg. 12711 for the parameters of the "universal" 
HMXB XLF, are shown for different values of SFR in the 
left panel of Fig|5| The right panel in Fig|^ shows the de- 
pendence of the most probable value of the luminosity of 
the brightest HMXB on the SFR of the host galaxy, de- 
scribed by ea. llAlll . along with its intrinsic 67% uncertainty. 
Filled symbols in Fig|^are the luminosities of the brightest 
so urce observed in sta r forming galaxies from the sample 
of lCrimm et alj i2003L and references therein). As is clear 
from the plot, the large difference in the maximum luminos- 
ity between low- and high-SFR galaxies, e.g. between the 
Milky Way and the Antennae galaxies, can be naturally un- 
derstood in terms of the properties of the probability distri- 
bution p (Z/max). No additional physical processes, affecting 
HMXBs formation in starburst galaxies, need to be invoked. 

4.3 Intermediate mass black holes 

The hypothetical intermediate mass black holes, probably 
reaching masses of ~ 10^~^Mq, might be produced, e.g. via 
black hole merges in dense stellar clusters, and can be asso- 
ciated with extremely high star formation rates. To accrete 
efficiently they should form close binary systems with nor- 
mal stars or be located in dense molecular clouds. It is natu- 
ral to expect, that such objects are significantly less frequent 
than ~stellar mass black holes. The transition from ~stellar 
mass BH HMXB to intermediate mass BHs should mani- 
fest itself as a step in the luminosity distribution of com- 
pact sources fFig llOl left panel) . If the cut-off in the HMXB 
XLF, observed at log(Lcut) ~ 40.5 corresponds to the max- 
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Figure 10. Illustration of the effect of hypothetical intermediate mass black holes on the Ljf — SFR relation. Left: The luminosity function 
of compact sources at different levels of star formation rate. Right: Corresponding Lx — SFR relation. The thin straight line shows the 
linear dependence. 



imum possible luminosity of "ordinary" ~stellar mass black 
holes and if at L > Lcut a population of hypothetical in- 
termediate mass BHs emerges, it should lead to a drastic 
change in the slope of the Lx-SFR relation at extreme val- 
ues of SFR (Fig llOl right panel). Therefore, observations of 
Lx— SFR relation for distant star forming galaxies with very 
high SFR might be an easy way to probe the population of 
intermediate mass black holes. 



4.4 Low mass X-ray binaries 

4.4-1 Lx — stellar mass relation and maximum luminosity 

As was shown by iGilfanovl ll2004) the luminosity distribu- 
tion of low mass X-ray binaries in nearby early type galaxies 
and bulges of spiral galaxies can be described by a "univer- 
sal" XLF whose shape is approximately the same in different 
galaxies and whose normalization is proportional to the stel- 
lar mass of a galaxy. The shape of the "universal" LMXB 
XLF is significantly more complex than that of HMXBs. 
It appears to follow the power law at low luminosi- 
ties, gradually steepens at log(Z/x)^37.0 — 37.5, and has 
a rather abrupt cut-off at log{Lx) ~ 39.0 — 39.5. In the 
log(I/x) ~ 37.5 — 38.7 luminosity range, it approximately 
follows a power law with the differential slope of ~ 1.8 — 1.9. 

Given the shape of the XLF, the total luminosity of 
LMXB sources in a galaxy is defined by the sources with 
log(I/x) ~ 37 — 38, the contribution of the brighter and, 
especially, weaker sources being less significant. Therefore 
the non-linear regime in the Lx — M, relation, although 
does exist for log(M*) < 10.0 — 10.5, is less important than 
in the Lx— SFR relati on for high ma ss X-ray binaries (see, 
for example, Fig. 14 in lGilfanovll2004il . 

Significantly more pronounced is the dependence of the 
luminosity of the brightest source on the LF normalization, 
i.e. on the stellar mass of the host galaxy. In order to study 



this dependence we used the broken power law approxima- 
tion for the LMXB XLF from lGilfanwl (|2004) and performed 
a series of the Monte-Carlo simulations, similar to those de- 
scribed in section|21 The probability distribution of the max- 
imum luminosity p (Lmax), obtained from these simulations, 
is shown in the left panel in Fig llll for different values of the 
stellar mass of the host galaxy. The right panel shows the 
dependence of the most probable value of the maximum lu- 
minosity and of its 67% intrinsic spread on the stellar mass. 
Its broken line shape is a result of the broken power law 
approximation of the LMXB XLF used in the simulations. 
Solid symbols show the observed values of the maximum lu- 
minosity for the number of nearby early type galaxies, bulges 
of spiral g alaxies and for L MXBs in the Milky Way from the 
sample of IGilfanovl (|2004| . and references therein). 

Similarly to HMXBs in star forming galaxies it is obvi- 
ous from Fig llll that the significant difference in the value 
of the luminosity of the brightest source can be naturally 
explained by the properties of the luminosity function of 
LMXBs. The same effect leads to an artificial (unphysical) 
dependence of the average luminosity of low mass X-ray bi- 
narie s in a galaxy on its stellar mass (e.g. Fig. 17 in IGilfanovl 
l2004) . So far there is no evidence for a significant change of 
intrinsic properties of low mass X-ray binaries with galactic 
environment. The difference between the luminosity of the 
brightest LMXB in massive elliptical galaxies and the bulges 
of spiral galaxies can be understood based on the probability 
arguments. 

4-4-2 Variability 

As in the case of HMXBs, in the linear large-mass limit, 
log(Af,) > 10.5, the fractional rms of the aperiodic vari- 
ability of the combined emission of LMXBs follows the 
rms tot oc 1 / V * averaging law. Owing to the shape of 
the LMXB XLF it decreases rather quickly with the stellar 
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Figure 11. Left: Probability distribution for the luminosity of the brightest LMXB for different values of the stellar mass, computed 
from eg. 1271 . using the parameters of the broken power law approximation to the "universal" luminosity function of LMXBs. Right: 
Expected luminosity of the brightest LMXB source vs. stellar mass of the host galaxy. The thick solid line shows the most probable 
value of the Lmax , the shaded area shows its 67% intrinsic dispersion, both obtained from the probability distribution given by eq. 1271 . 
The filled circles show maximum observed luminosities of LMXB sources in the Milky Way and several nearby galaxies, studied with 
Chandra. The "broken" shape of the predicted dependence is a consequence of the broken power law approximation of the "universal" 
LMXB XLF, used in the calculations. 



mass of the galaxy in the non-linear low-mass regime as well 
fFig |12t . Consequently, in massive elliptical galaxies, with 
stellar mass log(M*) ~ 11.0 — 11.5, the fractional rms vari- 
ability of the total emission will be suppressed by a factor of 
~ 10—15 with respect to the rms of individual sources. In a 
galaxy similar to the Milky Way with log(M.) ~ 10.5-10.7 
the suppression factor is ~ 5. Considerable variability on 
the level of~l/4 — l/2of that of individual X-ray binaries 
can be expected only for light bulges of spiral galaxies with 
masses in the log(M,) ~ 9.5 — 10.5 range. 

Fie ll3l compares the dependence of the fractional rms 
on the most probable value of the total luminosity for high 
and low mass X-ray binaries. In the bright luminosity end, 
log(Ljs:) >39.5, the X-ray emission from early type galaxies 
is expected to be significantly, up to a factor of ~ 7, less 
variable than from star forming galaxies. 



example, Lj, are the luminosities of individual X-ray bina- 
ries in a galaxy and Ltot is its total X-ray luminosity due to 
HMXB population. In this example the normalization of the 
luminosity function, i.e. the number of HMXBs in the galaxy 
n, is proportional to its star formation rate. We showed that 
due to statistical properties of the probability distribution 
p{Ltot) the result of a measurement of the total luminos- 
ity of a randomly chosen galaxy might deviate significantly 
from the intuitively obvious expression 

(Ltot) = / L — dL(xn(x SFR (32) 

These properties of p(I/tot) can result in a surprising non- 
linear dependence of the total luminosity on n. They can 
also cause anomalous variability of the combined emission 
in an apparent violation of the rms oc l/\/n averaging law 
for uncorrelated variations of individual sources. 



5 SUMMARY 

We studied the statistical properties of the combined emis- 
sion of a population of discrete sources. Namely, we consid- 
ered the properties of their total luminosity 

Ltot =^Lk (31) 
fc=i 

and its dependence on the number of sources n or, equiva- 
lently, on the normalization of the luminosity function (LP). 
Using high mass X-ray binaries in star forming galaxies as an 



Our results can be summarized as follows: 

(i) The probability distribution p (Ltot) can be computed 
numerically from eq.JJOJ or alternatively ea. l|23|l . Examples 
are shown in FiglTI 

(ii) The relevant characteristics of the p(Ltot) distribu- 
tion are: (1) its mode Ltot - the value of the random vari- 
able Ltot for which p(Ltot) has the maximum and (2) the 
expectation mean (Ltot) defined by eg. 1321 . 

It is the mode of the p (Ltot) distribution that predicts the 
most probable value of the total luminosity of a randomly 
chosen galaxy. If many galaxies with ~same total number 
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Figure 12. Variability of the total emission of low mass X-ray 
binaries. Dependence of the ratio rms tot /rmso on the stellar 
mass of the galaxy. The thick solid line shows the most probable 
value of the rms tot /rmsQ, the shaded area shows its 67% intrinsic 
dispersion, both obtained from the Monte -Carlo simulatio ns for 
the "universal" luminosity function from iGilfanovl i2004) . The 
dashed line shows asymptotical behavior at large M, . The filled 
circles correspond to LMXBs in the Milky Way and several nearby 
galaxies, computed from eg. 1281 . using the observed luminosities 
of X-ray sources in these galaxies. 



of sources n are observed the measured values of their total 
X-ray luminosities are distributed according to the p(Ltot) 
distribution whose shape depends on the chosen value of n 
(Figini left panel). The average of the measured values of 
Ltot is equal to the expectation mean (itot) and is always 
proportional to n in agreement with eg. 1321 . 

(iii) For small values of the LF normalization, Ltot and 
(I/tot) do not equal each other (Fig0 and the Ltot — n re- 
lation is non- linear (FigEHSlEJ- Only in the limit of n > 1, 
Ltot = {Ltot) and Ltot depends on n linearly. The threshold 
value of n depends on the LF shape and can be arbitrarily 
large. 

(iv) The skewness of the p [Ltot) probability distribution 
in the non-linear regime (Fig0 results in an enhanced and 
significantly asymmetric intrinsic dispersion of the measured 
values of Ltot ('Fig l3l6ll . Its non-Gaussianity precludes the 
use of the standard fitting techniques in analyzing the Ltot — 
n relation (e.g. Lx~SFR relation for star forming galaxies), 
such as minimization technique. 

(v) The variability of the total emission (e.g. aperiodic 
variability of X-ray emission of a galaxy due to superposi- 
tion of variabilities of individual sources) in the non-linear 
regime decreases with the number of sources more slowly 
than rms cx l/^/n law for uncorrelated variations of indi- 
vidual sources, resulting in anomalously high variability of 
the total emission. In the linear regime the rms oc ^/^n 
dependence is restored (section |21 Figl^. 

(vi) The amplitude of the discussed effect depends on the 



Figure 13. Comparison of variability of low and high mass X- 
ray binaries. Dependence of the ratio rms tot /rmsQ on the total 
X-ray luminosity of the galaxy due to X-ray binaries. The thick 
solid lines show the most probable value of the rms tot /rmsQ, the 
shaded areas show its 67% intrinsic dispersion, both obtained 
from the Monte-Carlo simul ations for the respe ctive "univer- 
sal" lu minosity functions from lGrimm et al.l i20n,'J) and IGilfanovl 
J2004) . For both curves, the linear parts at high Lx follow 
rms oc l/v/Lx averaging law. 



shape of the luminosity function. For a power law LF it is 
strongest for 1 < a < 2, is unimportant for shallow lumi- 
nosity functions with a < 1, and gradually diminishes with 
increase of a above a = 2 fFig lTlOl . 

We illustrate these results using the example of the com- 
bined emission of X-ray binaries in galaxies and its depen- 
dence on the star formation rate and on the stellar mass of 
the host galaxy. 

(i) For the slope of the HMXB "universal" XLF, a ^ 1.6, 
the discussed effects are strongest with a significant non- 
linear regime in the Lx— SFR relation at SFR < 4— 5 A/0 /yr. 
The predicted Ljc — SFR dependence is in good agreement 
with observations (FigEJ. Given the shape of the "univer- 
sal" LMXB XLF no significant non-linearity of the Lx — Af« 
relations is expected, also in a good agreement with obser- 
vations. 

(ii) The Lx— SFR relation can be used to constrain the 
XLF parameters of HMXBs in distant unresolved galaxies 
including the galaxies observed with Chandra in HDF-N at 
redshifts z ~ 0.2 - 1.3 (sections gTHl EH FigslTIITnil. 

(iii) Both for high and low mass X-ray binaries a strong 
dependence of the luminosity of the brightest source on 
the SFR and stellar mass of the host galaxy is expected. 
The Lniax— SFR and Lmax — Af« dependences predicted from 
the respective "universal" XLFs explain well the results of 
Chandra observations of nearby galaxies (Fig|5] and lllll . 
The significant difference in the luminosity of the brightest 
LMXB between bulges of spiral galaxies and giant ellipticals 
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or between the brightest HMXB in the Milky Way and in 
starburst galaxies can be understood based solely on prob- 
ability arguments. 

(iv) We predict enhanced variability of X-ray emission 
from star forming galaxies due to HMXBs, significantly 
above the oc l/\Ai averaging law. For SFR<5 M0/yr the 
expected fractional rms of variability of the combined emis- 
sion of HMXBs does not depend on the star formation rate 
and approximately equals ~ 1/3 — 1/2 of the rms of in- 
dividual sources (Fig|5J. On the contrary, variability of X- 
ray emission from early type galaxies due to LMXBs will 
be significantly suppressed because of the averaging effect 
(Figini, upto ~ 3 - 10 times in the ~ 10^" - 10^^ M© 
stellar mass range. For the same total luminosity star form- 
ing galaxies are expected to have significantly larger frac- 
tional rms, than massive elliptical and SO galaxies assuming 
that fractional the rms of individual sources are comparable 

(FigESJ. 
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APPENDIX A: APPROXIMATE SOLUTION 
FOR THE TOTAL LUMINOSITY 

We consider the case of a power law LF, eg. 11251 . with a 
slope a > 0. The probability distribution for the luminosity 
of the brightest source in the sample of n sources is defined 



by eq.lHTJ. The maximum of this distribution gives the most 
probable value of the luminosity of the brightest sources: 



-£max = min (l/max , L2) 
-'-'max \ T t ^ / -I \ 

-— =14-— — n-1 



(Al) 



Similarly, the probability distribution of the minimum 
luminosity in the sample is: 



P(imiii) = [Pl{L > Lmin)]" ^Pl(Z/min) n 



(A2) 



Contrary to p (Lmax) p (imin) declines steeply at L > _Li for 
any n (for a > 0). Within the accuracy of this approximation 
we can assume p(Lmin) = 5{Luiin — Li), i.e. Lmin = Li. 

The total luminosity of n sources distributed between 
Li and Lmax according to the power law with the slope of 
a (a 7^ 1, a 7^ 2) can be approximated as: 



1 ^ r2 — a r2 — a 

Ltot ~ n 



2-a L 



1 — a r 1 — a 
max 



(A3) 



Knowing the probability distribution p (Lmax) the prob- 
ability distribution pn(Ltot) can be calculated as: 



(itot) ~ P 



(A4) 



where Lmax ~ Lmax (Ltot) is iuverse function to ea. ljA3|l and 
p(Lmax) is given by eq.JJlJ. 

The most probable value of Ltot is defined by the con- 
dition 

dpJ^ = (A5) 
"Ltot 

With eg. (1271 . IIA31 and 1A4II the above equation can be 
transformed to: ^ 

(q - 2)^''" - (q - 2)(1 + Q - n + «n)C'+" + (A6) 
+ {a - lf{l + n)C + [l + {a- l)n] = 



Li 



(A7) 



or, equivalently: 

^ (g - 2)e'° -K2 + « - Q^)g^+" + (q - i)'r - e 

(a-l)[(a-2)e+"-(a-l)e°+e] 
Q 7^ 1, Q 7^ 2. Because of the simplifying assumption 
Pmin(Lmin) — 5(Lmin — Li) and the approximate nature of 
eo. <A31 . the probability distribution defined by ea. llA41 . is 
valid only for Ltot < (Ltot) and is undefined otherwise. This, 
however, is sufficient for our purpose as Lmax ~ L2 corre- 
sponds to the break in the Ltot — n relation (Fig|5|and|^, 
above which Ltot = (Ltot). 



Al The practical recipe 

The Ltot — n relation can be computed parametrically using 
eos. llA71 and IIA31 The practical recipe is for a set of values 
of Lmax, Li < Lmax < L2 , to computc Tl from eg. 1 A 71 and 
Ltot from eg. lA31 . The pairs of values (Ltot,Ti) define the 



Note that in eq. IX6l and |^ Lmax is a parameter rather than 
the most probable value of the maximum luminosity. The latter 
is 



defined by the eg. lAll . which is exact. 
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Ltot — n relation before and up to the break. Above the 
break, Ltot = {Ltot) and can be computed from ea. lA3ll 
with Lmax = L2 and n > ribroak - free parameter. The n 
in the obtained Ltot — n relation can be transformed to the 
normalization A via eg. llH . 

The approximation defined by the eas. llA7t and JA3t is 
compared with the results of the exact calculation in Fig|21 
It is accurate within ~ several per cent everywhere, except 
the break region, where its accuracy is ~ 10 — 20%. 



A2 Asymptotics 

Using the approximate solution for Ltot from AppendixIXIwe 
consider the asymptotical behavior of the Ltot — n relation 
in the limit of L2/L1 — > cx, a > 1. 

From ea. llA7ll variable ^ is related to the number of 
sources by (see footnote |5J: 



(q- l)n + 0(l) 



Ln 



(A8) 



in the limit of n >> 1 or equivalently Lmax >> Li. Al- 
though in the limit n ^ 00 this approximation is valid for 
any a > 1, its accuracy deteriorates considerably for q < 1.6 
where it can be improved by replacing n in eg. IAS! : 



1 + Q - 



1 < Q<1.6 



(A9) 



(a-l)(«-2)' 
The most probable value of the total luminosity is given by 

Ltot _ l-[(Q-l)n]^ 



(Ltot) 1 - (La/Li)^-'^ 



(AlO) 



a > 1, a 7^ 2, n >> 1. Using ea. llllll . 1121 it can be ex- 
pressed via the normalization of the luminosity function A 
or transformed to the relations for Ltot- As with ea. HASII . 
the accuracy of eo. lAlOl l can be significantly improved us- 
ing eg. lAOt for a < 1.6. 

For 1 < Q < 2 the Ltot — n relation shows a sharp break 
between the non- linear and linear regimes (Fig |2l3ll . From 
eg. lAlOj one can obtain: 



Ltot oc 



n < ribrcak 
n > ribrcak 



(All) 



This is in an agreement with eq.® based on simple quali- 
tative arguments. The position of the break in the Ltot — n 
relation can be obtained from eo. lA7t substituting £, ~ 
and using the fact that ^ >> 1: 



^brcak 



(A12) 



Expressed in terms of the the normalization A of the lumi- 
nosity function eo. H25^ it is: 



Ah 



(A13) 



As intuitively expected, the break position expressed in 
terms of the normalization of the luminosity function does 
not depend on the low luminosity cut-off Li and is defined 
only by the slope of the luminosity function and the high 
luminosity cut-off (see discussion in section [2.2.211 . The to- 
tal luminosity at the break, however, depends on the low 



luminosity cut-off for steep luminosity function with a > 2: 

if 1 < a < 2 



_t2_ 

2-Q 



^x(^)"-^ ifa>2 



(AM) 



as the total luminosity for a > 2 is defined by the sources 
near the low luminosity cut-off. 



APPENDIX B: VARIABILITY OF THE TOTAL 
EMISSION 

In the linear regime of Ltot — n relation, the fractional rms'^ 
of the collective emission is inversely proportional to the 
number of sources: 



(1 — a)(3 — a) n 
(a-l)(3-a) \Li) 



if a < 1 
i if 1 < a < 2 



(2^1^ (J^y-^l if2<a<3 

(a — 1)(3 — ck) \ Li / n 



(Bl) 



(cK — l)(cK — 3) n 

or, equivalently, to their total luminosity 



if a > 3 




if a < 2 

if 2 < Q < 3 (B2) 
if a > 3 



The above formulae are valid in the limit Li << L2. 

Similarly to the most probable value of the total lumi- 
nosity (Appendix^ the fractional rms of the total emission 
can be approximately calculated substituting L2 in eg. 1291 
with some value Lmax < L2. In principle, the probability 
distribution for the rms tot in eg. 12811 could be derived using 
the probability distribution for Lmax- The maximum of this 
probability distribution would give a sufficiently accurate 
approximation for rms tot /rmso- However, for simplicity we 
use the value of Lmax from eg. lA7ll - In the limit L2 >> Li 
one finds for 1 < a < 3: 



(a -2)2 



rmsl (a- l)(3-a) (^2-° - 1)2 n ^^^^ 

for n < ribroak, where ^ is defined by eg. lASII with substi- 
tution of ea. lA9ll for 1 < a<1.6 and ribrcak is defined by 
eg. lA12L For a < 1 and a > 3 there is no non- linear regime 
and the fractional rmstot obeys eg. lBH and IB2II for any n. 
From eq.lESJ one finds for 1 < a < 2 



(3 -a) 



n < ribroak 



(B4) 



i.e. in the non-linear regime the fractional rms of the collec- 
tive emission does not depend on the number of sources. 

The accuracy of egs. lB3ll and 1B4I I is sufficiently good 
for a > 1.5 but deteriorates for smaller values of a. In the lin- 
ear regime eqs. dHTl and (ISU are almost precise, their only 
approximation is in neglecting higher orders of the L1/L2. 



